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Abstract 

Geometric  paring  (GP)systems  are  being  developed  by  AOTC  Ft.  Hood  to 
augment  or  replace  laser-pairing  instrumentation  for  long-range  and  unconventional 
weapon  tests.  Such  systems  require  terrain- augmentation  to  rapidly  calculate  accurate 
line-of-sight  (FOS)  algorithm  accommodating  vegetation  and  obscurants  and  non  linear 
trajectory  prediction  through  complex  near  earth  terrain. 

At  one  meter  resolution,  both  vegetation  types  and  terrain  elevation  uncertainty 
complicate  the  FOS  calculations.  We  will  present  advances  in  statistical  FOS  algorithm 
for  calculating  LOS  probabilities,  which  accounts  for  elevation  measurement  errors  and  a 
varying  vegetation  density  along  the  LOS  path.  The  use  of  PVNT  deterministic  partial 
visibility  calculations  is  presented  as  low  cost  algorithm  test  methodology.  Then  the 
results  of  several  accuracy  vs.  speed  studies  using  this  method  are  presented. 

In  addition  requirement  for  small,  rugged,  lightweight  equipment  capable  of 
being  unobtrusively  carried  by  dismounted  infantry  has  generated  the  necessity  to 
develop  wide-area,  1 -meter  resolution,  compact  terrain  database  formats  and  LOS 
algorithms  suitable  for  execution  in  small  commercially  available  hand  held  devices.  We 
will  describe  a  Heirarchical  LOS  culling  technique  designed  to  increase  the  speed  of  LOS 
calculations  suitable  for  execution  in  small  devices. 

Keywords:  line-of-sight,  geometric  pairing,  terrain  database,  battlefield  digitization, 
elevation  uncertainty 
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1.  Introduction 


The  Mobile  Automated  Instrumentation  Suit  (MAIS)1  is  an  advanced 
instrumentation  system  that  supports  the  testing  of  U.S.  Army  weapon  systems  and  tactical 
doctrine.  MAIS  provides  a  realistic  battlefield  environment  for  simulated  force-on- force 
engagements,  real-time  casualty  assessment  (RTCA),  and  reconstruction  of  events  for 
posttest  analysis.  RTCA  requires  knowledge  of  who  shoots  at  whom,  that  is,  pairing. 

Substituting  laser  pulses  for  bullets  and  equipping  engagement  participants  with 

pulse  detectors  performed 
traditional  pairing,  as 
shown  in  figure  1 .  Though 
laser-pairing  systems  have 
served  the  operational  test 
community  well  for 
decades,  problems  in 
maintenance,  accuracy, 
safety,  mismatch  in 
obscurant  specific  bullet- 
versus-pulse  propagation 
characteristics,  and  its 
inability  to  characterize 
smart  and  guided  weapons 
systems  have  lead  to  the 
investigation  of 
alternatives. 

One  such  alternative 

is  the  Geo-Pairing  (GP)  system  being  developed  at  USAOTC  at  Ft.  Hood  Tx.  GP  systems 
use  global  positioning  system  (GPS)  position  and  pointing  sensors  as  input  in  order  to 
calculate  weapon-firing  trajectories.  .  Their  potential  advantages  include  low 
maintenance,  no  radiation  hazard,  potentially  higher  accuracy,  and  the  possibility  of 
including  advanced 
capabilities  in  the 
pairing  algorithms. 

Examples  of 
capabilities  that 
could  be  included 
are  detection 
algorithms  and 
ordinance-specific 
obscurant 
calculations.  In 
short  GP  enables  the 
MAIS  to  overcome 
most  of  the 
limitations  of  laser 

pairing,  that  is,  *  ordinance  specific  obscurant  calculation 


Fig  2  -Geometric  Pairing 
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Fig.  1  -  Traditional  Laser  Pairing 
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engaging  a  target  through  smoke,  rain,  fog,  and  foliage  and  at  longer  ranges  than  are  safe 
with  a  laser. 

In  order  to  fully  realize  the  potential  of  a  GP  system,  knowledge  of  the  terrain 
must  be  included  in  the  calculations.  Without  terrain  knowledge,  a  trajectory  calculated 

only  from  a  shooter’s 

_Fig  3  Terrain  Augmented 
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position  and  pointing 
angles  would  penetrate 
through  all  obstacles  in 
the  path.  This  would 
make  the  system  useless. 
Terrain  knowledge  at  a 
level  of  fidelity 
commensurate  with  the 
weapons  system  effects 
is,  therefore,  required  to 
calculate  a  realistic 
trajectory  on  an  actual 
battlefield.  Such 
trajectories  apply  to  the 
path  of  the  projectiles  as 
well  as  to  the  path  of 
photons  used  in  sensing 

devices  to  facilitate  detection,  aim,  and  fire  operations.  The  latter  is  referred  to  as  the 
LOS  problem,  and  literature  describing  the  Department  of  Defense  standards  and 
algorithms  are  available.3, 4 

The  theory  of  LOS  calculations  in  terrain  when  terrain  elevation  and  obscuring 
material  is  taken  into  account  was  presented  at  the  2003  ITEA  Symposium,  in,  Lihue, 
Kauai,  Hi.  Methods  were  presented  for  calculating  LOS  probabilities  when,  in  addition  to 
elevation,  the  terrain  is  characterized  by  vegetation  height  (vgh),  vegetation  identification 

(vid),  and  an 

Fig  4  PVNT  Source  Databases 
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the  lsq.  m  element  encountered  on  the  top  of  the  ground  material  or  bottom  of  the  sky. 
Databases  containing  the  four  parameters  have  so  far  bee  been  built  in  five  locations 
shown  on  the  map  insert  of  Figure  4.  These  database  were  built  by  the  Perspective  View 
Nascent  technologies  (PVNT)  described  in  the  sister  paper  presented  at  this  conference. 
PVNT  is  a  low  cost  PC  based  systems  designed  to  address  the  generation  and  utilization 
of  1  meter  resolution  terrain  in  large  area  tactical  battlefield  simulations.  The  Naval 
Postgraduate  School,  TRADOC,  OTC  Ft.  Hood  and  Nascent  Systems  Inc  have  developed 
the  package  as  a  research  project.  It’s  major  advantages  are  that  it  runs  on  a  Laptop  PC, 
contains  an  extremely  efficient  ray  trace  algorithm,  performs  real  time  rendering,  handles 
complex  surface  features  for  low  angle  oblique  views,  and  is  freely  available  to  DOD 
users. 

It  should  be  pointed  out  that  even  with  tools  available  and/or  being  developed  the 
construction  of  metrically  accurate  databases  at  the  one-meter  resolution  level  over  such 
large  geographic  areas  is  an  expensive  undertaking  that  has  generally  been  funded  only 
for  the  parameters  actually  required  for  specific  operational  tests.  Hence  not  all 
parameters  are  accurately  available  for  all  areas.  Inclusion  of  vegetation  layer  and 
underlayment  penetration  probabilities7  required  for  probabilistic  LOS  are  ongoing 
efforts  and  parameters  for  additional  regions  will  be  recorded  as  funding  constraints 
permits. 

The  Probabilistic  LOS  algorithms  have  been  implemented  in  a  PC  104  based 
terrain  server  and  tested  at  Ft.  Hood  when  deployed  on  vehicles  and  tanks.  The  timing 
constraints  of  the  MAIS  system  requires  that  an  LOS  calculation  be  performed  in  less 
than  10  millisecond  including  request  message  and  response  message  overhead. 
Additionally  the  required  geographic  areas  over  which  such  systems  are  expected  to 
operate  have  been  expanded  to  include  air-to-ground  and  long  range  smart/guided 
weapons  systems.  On  top  of  that  the  host  system  in  which  such  LOS  algorithms  operate 
should  be  reduced  in  size  in  order  to  allow  dismounted  infantry  to  be  instrumented 
without  excessive  weight  burdens.  This  has  lead  to  the  development  of  new  LOS 
algorithms,  which  utilize  smaller  memories,  execute  faster,  and  address  larger  geographic 
areas  than  previously  possible.  The  remainder  of  this  paper  presents  the  design  and 
implementation  details  of  these  new  LOS  algorithms  capable  of  supporting  the  Terrain 
Augmaneted  Geometric  Pairing  mission. 


2.  LOS  Algorithm  in  Uncertain  Terrain 


The  theory  and  mathematical  formulation  of  the  basic  LOS  algorithm  capable  of 
handling  multiple  vegetation  layers  in  uncertain  terrain  was  presented5  at  the 
Test&Evaluation  of  Advanced  Technology  Systems  Conference  ,  Nov.  17-21,  2003, 
Lihue,  Kauai,  Hi.  This  section  will  summarize  the  algorithm  and  present  some  test  results 
methodologies.  Presented  results  will  provide  the  rational  for  a  next  generation 
hierarchical  design  intended  to  minimize  the  execution  speed  requirements  and  allow  the 
algorithm  to  run  over  large  areas  in  order  to  make  it  suitable  for  long  range  engagements 
encountered  in  modem  weapons  testing. 

2.1)  Basic  LOS  Probability  Algorithm: 

The  geometry  of  the  situation  is  exemplified  in  figure  5.  Here,  the  terrain 
correlation  cell  is  divided  into  smaller  steps  of  size  As,  and  the  ray  is  slanted.  At  each 
step,  the  composition  of  the  terrain  layers  changes.  Hence,  the  penetration  probability 
through  the  cell  is  no  longer  horizontally  homogeneous.  If  there  were  no  terrain  error, 
the  ray  would  enter  the  correlation  cell  at  101  meters,  travel  through  two  steps  of 
vegetation,  and  hit  solid  ground  in  the  third  step  at  elevation  99  meters. 

The  penetration  probability,  Pp(t  ,s),  of  penetrating  to  step  s  from  the  observer  in 
the  correlation  cell  t  is  the  product  of  the  individual  penetration  probability  through  each 
step. 
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The  most  useful  form  for  our  algorithm  is  the  recursive  form;  hence,  the  final  expression 
of  the  LOS  probability  is  given  in  terms  of  calculation  step  size  as: 


Eq.  1-  Losp(s)  =  Losp((int)  s  *  As  /  At) 

+oo 

*  fp(z(,), -.m,  01(1))  *2IjS/“W  " <’> ■  zr(she>  >de 

-oo 


where:  e  =  the  terrain  error  equal  to  z  -  zt 

s  =  length  in  step  sizes  from  viewer  to  target 
As  =  step  size  inside  the  correlation  cell 
t  =  correlation  cell  distance  from  viewer  to  target 
At  =  size  of  a  terrain  correlation  cell 
Cos(el)  =  the  cosine  of  the  elevation  angle 
zt(t)  =  elevation  of  terrain  as  a  function  of  terrain  correlation  cells 
at(t)=  error  of  terrain  as  a  function  of  terrain  correlation  cells 
Sm(e(s),  n(s),  zr(s)  =  the  half  penetration  distance  as  a  function  of  easting  e(s), 
northing  n(s),  and  ray  altitude  zr(s) 

£  =  the  sum  to  be  carried  out  from  the  first  step  in  the  correlation  cell 
s  =  t  *  At  /As  to  the  step  s  inside  the  correlation  cell 

2.1.1  Algorithmic  Explanation 


To  help  visualize  the  algorithm  being  described  lets  look  again  at  the  diagram  of 
the  basic  LOS  problem  shown  in  figure  5.  The  center  dark  line  is  the  nominal  ground 
value.  The  colored  bands  are  error  bound  bins  each  of  which  corresponds  to  a  probability 


The  algorithm  has  calculated  the  probability  of  the  ray  reaching  the  left  side  of  this 
terrain  correlation  bin.  Lets  say  this  is  LOS_T  =  .75. 

In  the  first  step  the  ray  is  three  bins  high.  In  order  for  the  ray  to  be  stopped  the 
terrain  must  have  a  +  3  bin  error.  If  this  were  the  case  the  ray  would  actually  go  through  a 
vegetation  layer  with  a  step  penetration  probability  of  .5.  The  mask  shown  as  a 
transparent  rectangle  shows  the  probability  values  for  penetrating  this  step.  All  are  1 
except  the  top  hence  the  probability  of  penetrating  this  far  is  the  mask  multiplied  by  the 
terrain  error  probability  added  together.  This  is. 5*. 02+  09+.22+.29+.22+.09  =  .92.  We 
then  multiply  this  by  the  probability  of  entering  the  terrain  correlation  area  to  get  the  line 
of  sight  probability  for  penetrating  all  the  way  through  step  1 . 

Next  the  second  step  is  taken.  The  ray  is  lower  and  the  left  side  rectangle  shows 
the  step  penetration  mask.  Since  the  ray  is  lower  the  terrain  errors  above  +1  bin  would 
make  the  ray  invisible  hence  the  two  upper  bins  get  zeros.  The  cumulative  probability 
(PC)  of  penetrating  through  the  second  step  is  shown  on  the  right  rectangle  and  is  simply 
a  multiplication  of  the  current  penetration  probability  with  the  previous  cumulative 
probability. 

The  third  step  the  ray  is  one  bin  above  the  nominal  vegetation  layer.  If  a  terrain 
error  of  one  bin  is  made  the  penetration  probability  is  .5.  If  the  error  is  higher  the 
probability  is  0.  Hence  the  cumulative  probability  lowers  the  .5  by  one  bin. 

The  fourth  step  the  ray  is  still  in  the  same  bin.  Now  the  penetration  probability  of 
.5  if  a  plus  one  bin  error  is  made  is  multiplied  with  the  previous  cumulative  probability 
to  get  a  .25  for  this  bin. 

The  last  bin  the  ray  has  moved  down  to  the  vegetation  layer.  Hence  the  step 
penetration  probability  mask  has  moved  the  .5  down  to  this  bin.  If  no  errors  are  made  in 
the  nominal  terrain  altitude  then  the  ray  will  have  a  50%  chance  of  getting  through  this 
step.  The  cumulative  probability  of  getting  through  is  zero  unless  either  no  error  or  a 
negative  error  is  made  in  the  terrain.  The  probability  of  getting  through  the  five  steps  is 
.5*.29+.22+.09  =  .455. 

The  ray  is  partial  blocked  in  the  last  step.  If  the  terrain  had  fallen  here  then  the  Pp 
mask  would  have  ones  in  the  higher  bins  and  the  ray  penetration  probability  would 
remain  identical  to  the  previous  bin.  It  can  be  seen  that  the  penetration  probability  is 
determined  by  the  highest  obstruction  the  ray  passes  over. 

Probabilities  between  terrain  correlation  areas  are  AND’d  hence  if  LOS  T  was 
.75  to  get  to  the  left  side  of  the  probability  of  going  through  all  the  way  to  the  right  is 
.455*. 75  =  .34125. 

2.2)  Algorithm  Characteristics 

Implementing  Basic  LOS  Probability  Algorithm  requires  the  definition  of  the 
variables  defined  in  table  2-1.  Shown  are  the  program  variables  in  the  first  column,  the 
definition  in  the  second  column,  and  the  mathematical  symbol  from  equation  1  in  the  last 
column.  We  are  assuming  a  four  layer  terrain  includes  sky,  vegetation,  under  cut ,  and 
bare  earth.  The  height  of  these  layers  are  defined  in  the  database  for  each  1  meter  element 
the  penetration  value  applied  are  input  variables  set  at  runtime. 


S_HALF  Vegetation  penetration  half  distance(m)  Sm(e(s),  n(s),  zr(s) 

U_HALF  Vegetation  undercut  half  distance(m)  S\/2(e(s),  n(s),  zr(s) 

T_SIGMA  Terrain  sigma  error  value  at(t) 

DELTAT  Terrain  correlation  cell  length(m)  At 

DELTA_Z  Size  of  terrain  elevation  bins  ,  vertical  resolution  Not  def. 
PROBABILISTIC-EXIT  -  Stop  calculation  when 

Probability  is  less  than  this  value  Not  def 

Table  2-1  Main  LOS  Probability  Algorithm  Parameters 


Though  many  tests  were  performed  to  characterize  the  algorithms  with  changes  in 
input  values  the  two  most  significant  are  the  terrain  sigma  error  value  and  the  vertical 
resolution  of  the  probability  function  (see  P(z(t),  zt(t),  ot(t))in  figure  5).  The  significance 
of  the  terrain  error  can  best  demonstrated  using  a  high-resolution  terrain  simulator  PVNT, 
as  a  low  cost  alternative  to  field  trials  which  should  be  of  special  interest  to  the  test 
community. 


2.2)  PVNT  Partial  Visibility  Vs.  Probabilistic  LOS  Test 

The  PVNT  program  generates  line  of  sight  values  based  upon  a  deterministic  ray 
trace  algorithm.  The  default  places  a  2x2meter  rectangle  centered  at  the  target  point  and 
calculated  the  fractional  visibility  of  the  rectangle  from  the  view  point.  The  fractional 
visibility  is  defined  as  the  fraction  of  the  rectangle  visible  divided  by  the  total  area.  In 
practice  the  ration  of  visible  to  total  pixels  in  the  view  is  used. 


The  fractional  visibility  for  a  set  of  reference  test  points  calculate  in  PVNT  are  listed  in 


table  3-1  below.  Also  shown  are  the 
perspective  view  to  the  target 
rectangle  from  the  view  point.  Four 
test  points  are  from  a  helicopter 
position  at  point_21  (elevation 
362m)  while  two  more  are  from 
helicopter  position  at 
point_22(elevation  321m).  The 
points  were  chosen  because  the 
target  rectangles  were  close  to 
obscuration  and  hence  represent 
border  line  visibility  cases  which 
are  of  most  interest.  Easting  and  northing 
coordinates  are  shown  in  the  screen  print  out 
above  while  parameters  are  shown  in  the  box  on 
the  right. 

The  terrain  Sigma  value  was  varied  from  1 
meter  to  10  meters  to  show  the  effect  of  terrain 
error  on  probability  calculations.  In  order  to  get  a 
better  feel  for  the  meaning  of  the  probabilities  we 
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Table  2-2  PVNT  Fractional  Visibility  Vs.  Probabilistic  LOS  for  Various  Terrain 

error  sigma  values 


placed  these  values  along  side  the  perspective  views  of  the  target  rectangles  as  seen  from 
the  viewpoints.  These  are  shown  in  table  2-2.  The  right  column  shows  the  fractional 
visibility  and  the  perspective.  At  the  center  of  each  view  is  a  red  target  rectangle.  This 
can  be  partially  obscured  by  intervening  terrain. 

Comparing  the  fractional  with  the  probabilistic  calculations  one  can  notice  the 
expected  trends.  At  low  terrain  sigma  the  visible  points  in  general  are  very  visible  with 
probabilities  close  to  one.  As  the  terrain  becomes  uncertain  the  probable  visibility 
decreases  since  there  is  a  high  likelihood  that  some  terrain  point  is  high  enough  to  be  in 
the  way. 

For  grazing  point  cases  such  as  point_16  the  fraction  of  the  target  visible  is  23%  if 
the  terrain  is  perfect  as  modeled  in  the  database.  The  center  point  of  the  target  is  below 
the  intervening  trees  and  is  not  visible.  Hence  at  1  sigma  error  it  is  very  likely  that  the 
center  point  is  obscured  and  the  calculated  probability  of  visibility  is  11%.  As  the  terrain 
becomes  more  uncertain  the  probability  increases  that  the  center  point  is  visible.  Hence  at 
3  sigma  the  we  get  18%  closely  matching  the  fractional  visibility.  At  10  sigma  however 
the  probability  that  some  intervening  terrain  feature  causes  a  blockage  rises  and  the  LOS 
probability  moves  toward  zero. 

The  trend  for  probabilistic  LOS  to  go  to  zero  when  terrain  sigma  gets  large  is  due 
to  the  statistical  mathematics.  Some  terrain  element  will  statistically  stick  up  and  block 
the  view.  Though  correct  algorithmically  it  is  counterintuitive.  As  the  terrain  becomes 
uncertain  one  might  expect  the  LOS  probability  to  approach  50%. 

2.3  Speed  Vs.  Probabilistic  Accuracy 

A  four  kilometer  calculation  with  the  main  parameters  shown  in  the  last  section 
consumed  approximately  6  milliseconds.  The  distance  and  the  size  of  terrain  elevation 
bins  or  vertical  resolution  probability  parameter  DELTA_Z  .  In  the  next  section  we  will 
introduce  a  hierarchical  algorithm  to  reduce  the  run  time  dependence  on  distance.  Here 
we  note  that  the  vertical  resolution  probability  parameter  can  also  be  reduced  to  increase 
speed  at  the  cost  of  accuracy.  DELTA-Z  sets  the  elevation  integration  increment  for  the 
probability  numerical  integration.  Since  this  integration  is  performed  in  the  inner  loop 
decreasing  the  step  size  and  reduces  the  run  time  by  a  proportional  amount.  The  easiest 
way  to  trade  off  speed  with  accuracy  is  to  change  DELTA_Z. 

In  order  to  get  a  good  ground  trace  we  need  to  select  a  ray  path  which  is  long 
enough  to  show  results  and  along  which  the  ray  travels  in  a  grazing  direction  so  that 
reasonable  values  are  observable.  Values  selected  are 
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Between  these  two  points  the  ground  rises  gently  from  313  meters  at  the  view  point  to 
343  meters  at  the  target  location.  The  target  was  placed  3  meters  below  the  ground  in 
order  to  show  the  probability  calculation  as  the  ray  passes  through  the  ground.  Since 


there  are  errors  in  the  terrain  we  expect  the  probability  for  below  ground  rays  to  approach 
zero  but  not  be  zero  immediately  upon  striking  the  earth. 

The  following  table  2-3  shows  nine  calculations  performed  from  a  view  altitude 
320  meters  to  terrain  altitude  at  the  target  342  meters.  By  raising  the  target  from  2  meters 
below  the  ground  to  1  and  then  3  meters  above  the  ground  while  varying  DELTA  Z  the 
following  measurements  were  found. 


Relative  speed  ms 

Target  Altitude  -> 
DELTA  Z 

340 

343 

345 

18  ms 

.5 

.060 

.598 

.902 

9  ms 

1.0 

.080 

.664 

.933 

6  ms 

1.5 

.117 

.749 

.947 

Table  2-3  LOS  Probability  vs  Calculation  Speed 

The  execution  speeds  are  relative  and  do  not  imply  actual  numbers  over  the 
500meter  distance  but  rather  an  equivalent  4km  distance.  Deviations  in  the  worst  case  is 
.15  or  15%  of  the  probability  values.  This  is  an  extreme  case  encountered  only  for  very 
small  probability  calculations  when  the  ray  has  traveled  within  the  nominal  1  meter 
terrain  band.  Delivered  code  is  set  to  DELTA_Z=1.5.  We  believe  this  satisfies  the 
immediate  timing  requirements  for  MAIS  at  adequate  accuracies  however  clearly  at 
longer  ranges  will  prove  inadequate.  To  overcome  this  run  time  restriction  a  new 
hierarchical  LOS  algorithm  has  been  developed  as  described  in  the  next  section. 


3.  Hierarchical  LOS  Algorithm 

The  basic  LOS  algorithm  described  in  the  last  section  is  compute  intensive.  The 
longer  the  visibility  lines  between  the  viewer  and  the  target  the  longer  the  compute  time. 
Since  this  algorithm  must  execute  in  the  10msec  time  frame  in  order  to  fit  into  MAIS,  it 
is  important  to  develop  techniques  to  minimize  the  distance  dependence  on  execution 
time. 

In  order  to  make  the  algorithm  less  sensitive  to  view  distance  we  introduce  a 
multi-resolution  approach.  Simply  stated  the  LOS  calculation  can  be  done  using  64  meter 
step  sizes  in  order  to  provide  a  quick  and  approximate  answer  over  long  distances.  Such  a 
low  resolution  execution  will  identify  those  obscurant  points  along  the  path  around  which 
higher  resolution  calculations  need  to  be  performed. 


Figure  6  -  Hierarchical  LOS  Algorithm  Geometry 


Implementation  of  the  Hierarchical  LOS  Algorithm  requires  access  to  at  least  two 
terrain  resolution  databases.  In  our  case  the  low  resolution  database  consists  of  64  meter 
posts  while  the  high  resolution  post  size  is  either  1  or  2meters  on  a  side  depending  upon 
the  configuration.  The  low  resolution  only  requires  256  posts  per  sq.  km  and  large  areas 
of  up  to  64km  easily  fit  into  main  compute  memory  requiring  approximately  2Mbytes. 
However  low  resolution  data  will  require  several  Gigabytes  and  will  typically  not  fit  into 
small  machines  especially  hand  held  devices  intended  for  use  by  dismounted  infantry. 
Therefore  memory  management  for  large  areas  is  required. 

3.1)  Memory  Management  for  Large  Area  LOS 

The  geographic  LOS  database  problem  is  illustrated  in  figure  7.  A  pair  of  players 
labeled  To(T)  and  From(F)  are  separated  by  a  distance  along  which  the  LOS  is  to  be 
calculated.  In  general  terrain  data  along  this  distance  will  include  a  small  area  of  high 


resolution  data,  a  larger  area  of  low  resolution  data,  and  possibly  an  area  for  which  no 
data  is  available.  /p\ 


The  Hierarchical  LOS  Algorithm  avoids  distance  dependent  execution  times  by 
identifying  a  limited  number  of  high  resolution  blocks  in  which  the  calculation  is 
executed.  The  problem  is  that  we  do  not  A?)  know  which  high  resolution  blocks  are 
required.  If  there  were  no  time  constraint  /  the  obvious  solution  would  be  to  allow 

the  low-resolution 
calculation  to  identify  the 
high-resolution  blocks 
needed  and  retrieve  them 
from  disk  as  required.  Since 
it  takes  on  the  order  of 
milliseconds  to  retrieve  data 
from  disks  this  approach 
will  not  meet  the  timing 
requirements.  Hence  some 
form  of  preloading  is 
required. 

There  are  several 
algorithms,  which  may  be 
used.  The  simplest  is  to 

Figure  8  -  Multi  Resolution  LOS  Terrain  Geometry  ^oac^  a  rectanS^e  °f  high- 


High 

resolution 

Blocks 


resolution  data  around  the  target  player  location.  The  size  of  the  rectangle  can  expand 
depending  upon  the  amount  of  available  memory.  A  better  approach  is  to  load  only  those 
blocks  likely  to  contain  an  LOS  ray  obscurant.  Figure  8  shows  a  hypothetical  horizon 
contour  map  calculated  around  a  potential  target  T.  The  outer  line  is  the  horizon  as  seen 

from  T  while  the  inner  line 
surrounds  a  depression  into  which 
the  T  player  cannot  see.  A  small 
hilltop  is  also  visible.  Each  of 
these  visibility  contour  lines  can 
be  enclosed  in  a  high-resolution 
terrain  block.  All  the  blocks  taken 
together  then  represent  the  total 
space  within  which  LOS  ray 
blockage  is  likely  to  occur.  The 
actual  LOS  ray  calculation  from  a 
viewer  F  to  a  potential  target  T 
can  then  only  be  obscured  in 
those  areas  covered  by  the  loaded 
blocks. 

So  for  our  hypothetical 
example  a  64-meter  LOS 
calculation  would  identify  the 
hill(l),  the  ridgeline(2),  and  the 
area  immediately  around  the 
target  as  potential  areas  in  which 
high  resolution  calculation  should 
be  performed. 

A  block  diagram  for  the 
load  algorithm  is  shown  in  figure 
9.  The  flow  of  logic  assumes  the 
algorithm  is  run  in  a  thread  that 
waits  for  new  self  coordinates. 
Once  triggered  a  low-resolution 
view  shed  calculation  is  made  and 
the  number  of  blocks  from  the 
total  list  of  viewable  blocks  fitting 
into  main  memory  are  then  loaded 
from  the  disk.  Several  culling 
criteria  can  be  employed. 

Distance  from  the  target  is  the 
simplest  and  most  effective.  The 
high-resolution  blocks  are  then 
ready  for  an  LOS  calculation 
when  required. 


The  most  important  characteristic  of  the  memory  management  scheme  is  its  block 
structure.  Data  is  not  stored  in  large  arrays  but  in  blocks  covering  256x256  meters  in 
area.  Thus  a  map  is  required  to  translate  internal  block  location  into  easting  and  northing 
coordinates.  The  PVNT  terrain  database  are  already  in  the  required  block  structure  and 
resolutions  from  1  to  64  are  available.  Interfaces  are  therefore  built  directly  to  the 
existing  databases  as  shown  in  figure  10. 


4.  Conclusion 

The  Basic  LOS  Probability  Algorithm  allowing  calculations  in  uncertain  terrain 
with  multiple  vegetation  layers  requires  too  much  CPU  time  for  long  distance  pairing.  By 
introducing  a  multi  resolution  scheme  and  performing  LOS  ray  traces  at  the  64  meter 
block  size  this  problem  can  be  eliminated.  This  makes  high-resolution  1  meter  terrain 
calculations  feasible  over  large  playing  areas.  In  order  to  fit  the  required  large  amounts  of 
high  resolution  data  into  limited  memories  available  for  dismounted  applications  a  block 
structured  terrain  data  scheme  is  used  and  only  the  required  visible  areas  are  loaded.  This 
approach  addresses  both  the  high  speed  performance  for  large  playing  areas  in  limited 
memory  hardware  used  by  dismounted  players. 
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